Evidence for diquarks in lattice QCD 
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Diquarks may play an important role in hadron spectroscopy, baryon decays and color super- 
conductivity. We investigate the existence of diquark correlations in lattice QCD by considering 
systematically all the lowest energy diquark channels in a color gauge-invariant setup. We measure 
mass diflterences between the various channels and show that the positive parity scalar diquark is 
the lightest. Quark-quark correlations inside the diquark are clearly seen in this channel, and yield 
a diquark size of 0(1) fm. 
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Diquarks were originally proposed several decades 
ago 1] as a natural way to explain the absence of a large 
number of exotics otherwise predicted by QCD. Recently, 
it has been realized that they provide a natural explana- 
tion for an exotic baryon antidecuplet, the 0+ [3], that 
can not be accomodated in the quark model. The first 
evidence for the 6+ was reported by the LEPS Collabo- 
ration , but subsequent experiments cast doubts on its 
existence Q • Independently of whether a pentaquark ex- 
ists or not, diquark correlations may play an important 
role in the description of quark distributions and frag- 
mentation functions and in explaining the systematics 
of non-leptonic weak decays of light quark hadrons 
In addition diquarks are the central ingredient of cold, 
dense matter where they condense to form a color super- 
conductor. Despite their potential role in explaining such 
a variety of phenomena in hadronic physics, quantitative 
analyses that can directly assess whether diquark corre- 
lations are present in QCD are lacking. Lattice QCD is 
the method of choice for studying hadronic states. How- 
ever, diquarks are colored objects and need an appropri- 
ate formalism to study them using lattice simulations. 
One way is to fix the gauge as done in an early study 
on the mass of the diquarks [s^ . In this work we present 
a gauge invariant formalism, where we create color sin- 
glet states by considering diquarks in the background of 
a static quark This enables us to extract mass 

differences between baryons containing diquarks in the 
various channels. These mass differences, unlike diquark 
masses themselves, are well-defined and gauge invariant. 

One gluon exchange indicates quark-quark attraction 
in the color antitriplet 3c, flavor antisymmetric 3/ and 
spin singlet positive parity channel. Diquarks in this 
channel are referred to as "good" diquarks Q. The aim 
of this work is to check whether QCD dynamics supports 
attraction in this channel and compare it with other di- 
quark channels. Possible diquark configurations are cre- 
ated by two quark operators and insertions of the co- 
variant derivative. In this work we consider only diquark 
configurations with no derivatives since these are lower in 
energy. We therefore consider all 16 diquark multiplets 
that can be created by operators bilinear in the quark 



fields of the form q^CTq with C ~ «7o72 the charge con- 
jugation operator, and F = 1, 7^, 75, 757^, a^^. The 
positive parity channels are q^Cj^q and q^C^^^^q with 
spin zero and q'^Cjiq and q^Caoiq with spin one. The 
negative parity channels are q"^Cq and q'^Cjoq with spin 
zero and q^C'j5%q and q^Caijq with spin one. They 
create states that vanish in the non-relativistic limit and 
are excluded in quark models. We expect them to be 
higher in energy. 

In addition to mass differences, we probe correlations 
directly by examining the spatial distribution of the two 
quarks forming the diquark in the background of the 
static quark p. Correlations which persist as the static 
quark is removed are intrinsic to the diquark. We study 
these correlations using gauge invariant two-density cor- 
relators, used before 9] to probe hadron structure: 

Cr(r„,rd,i) = (0|Jr(O, 2t) (r„, t) J,?(rd, t) ^ (0, 0)|0) (1) 

where j/(r, t) /(r, i)7o/(r, t) :, / = it, d and 

Jr(x) = e"''^ u'^aix)Crdbix)±d^aix)CTu,,{x) s,ix) 

(2) 

where the -I- (-) sign corresponds to the flavor symmetric 
(antisymmetric) combination, and Sc denotes the static 
quark. Latin indices denote color. We use two degener- 
ate flavors of Wilson fermions on quenched configurations 
generated at /? = 5.8, 6.0 and 6.2 corresponding to lattice 
spacing a = 0.136, 0.093 and 0.068 fm as determined from 
the static quark force [l^l- At each lattice spacing, we 
perform measurements at three values of the pion mass 
in the range 570 — 910 MeV. By comparing results ob- 
tained at different lattice spacings but at the same pion 
mass we can assess discretization effects. In most cases 
the quenched approximation is good when the mass of 
the pion is higher than 600 MeV, but is expected to fail 
in the chiral limit. In order to check the validity of our 
results in this range of pion masses, we repeat our mass 
splitting measurements on a set of configurations with 
two degenerate flavors of dynamical Wilson fermions at 
/3 — 5.6 [lll|. In Table [1] we collect the parameters of our 
calculation. 
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TABLE I: Summary of our simulations, including the value 
of K, which determines the bare quark mass, the pion mass, 
m^, and nucleon mass, Afjv, in lattice and, for the latter, in 
physical units and the number of gauge configurations. 







aMN [Mn GeV] 


number of confs 


Quenched 16^ x 32 (3 = 5.8, 


a — 1.4/ Kjev 


0.1560 


0.619(2) 


1.139(8) [2.04(1)] 


OOD 


0.1575 


0.549(2) 


1.052(9) [1.55(1)] 


lOO 


0.1590 


0.473(2) 


0.961(9) [1.41(1)] 


zUU 


Quenched 16^ x 32 13 = 6.0, 


a-^ = 2.15 GeV 


0.153 


0.423(1) 


0.783(8)[1.68(2)] 


364 


0.154 


0.366(2) 


0.716(7)[1.54(2)] 


503 


0.155 


0.300(2) 


0.634(8)[1.36(2)] 


287 


Quenched 20^ x 40 13 = 6.2, 


= 2.94 GeV 


0.1510 


0.286(3) 


0.563(8) [1.66(2)] 


116 


0.1520 


0.214(3) 


0.487(9) [1.43(3)] 


166 


0.1523 


0.188(3) 


0.460(11)[1.35(3)] 


157 


Unquenched Wilson 24^ x 40 (3 = 


= 5.6, a-^ = 2.42 GeV 


0.1575|0.270(3) 


0.580(7)[1.40(2)] 


185 



The mass of our static-light-light baryon can be 
obtained from the large time limit of the correlator 
Gr{t) — {Jr{x,t)jl(x,0)). The effective mass, meg{t) = 
— \og[Gr{t)/Gr{t — 1)], becomes time independent for 
large t when the Ughtest state dominates (plateau re- 
gion), yielding its mass. To reduce the statistical errors, 
we need to keep t small by isolating the ground state 
as fast as possible. This is accomplished by suppressing 
higher excitations using Wuppertal smearing [13 ] on the 
source and sink with hypercubic (HYP) smeared spatial 
links [11] for the Wuppertal smearing function [ij]. In 
addition we use HYP smearing on the temporal gauge 
links U4{x) that enter in the construction of the static 
propagator given by 



Sstat.(x2,t2,Xi,ti) = 5 (X2-Xi) 



1 -I- 70 



n '^4(xi,t) 



(3) 

for t2 > ti. The exponential prefactor exp(— TOq(t2 — ^i)) 
has been dropped, since this introduces a constant shift 
in all energies by the bare heavy quark mass m^, which 
cancels in the mass differences that we measure between 
baryons with diquarks in different channels. 

We show in Fig. [T] the effective mass differences, Amoff, 
between the "good" (scalar) diquark created by taking 
r = 75, and the other channels. As expected the mass 
difference is zero for the other scalar diquark with posi- 
tive parity created by taking F = 7075- The other pos- 
itive parity channels, i.e. the "bad" (vector) diquarks, 
are degenerate and are clearly higher in energy than the 
"good" diquarks as pointed out in Ref. [8]. The nega- 
tive parity states, as expected, have even larger energies. 
Their effective masses are much noisier making it difficult 
to identify a clear plateau void of contamination from yet 
higher states. The general behavior is however similar to 
the positive parity diquarks: the vector channels are also 
degenerate and tend to have a larger mass than the scalar 
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FIG. 1: Effective mass difference, Arricff, in the various di- 
quark channels, a,t (3 = 6.0 for our lightest quarks. The ar- 
rows (from left to right) indicate the time slice of the source, 
of the density insertions and of the sink. 

channel created by taking F = /. To verify that the mass 
difference seen between the two scalar channels is signif- 
icant one would need to improve the quality of the mass 
plateaus. 

Effective color-spin Hamiltonian arguments [8j] predict 
that for heavy constituent quarks the mass difference 
Am = m-y; — rn-yj between the "good" and "bad" di- 
quarks scales like 1 /ruq-^ ruq^ , where the masses 
of the constituent quarks. On the other hand in the light 
quark regime we have cx m^. Therefore for intermedi- 
ate values of mq where both relations hold approximately 
we make the Ansatz Am cx In the opposite limit 
ruq 0, Am is expected to approach a constant value. 
These two behaviors can be connected by the Ansatz 

ci 



Am 



1 + (m^/c2)*' 



(4) 



with ci and C2 to be determined from the lattice data. 
We show in Fig. [2] results for Am as a function of m^ 
for our three /3 values. All lattice data should fall on 
a universal line if we are close to the continuum limit. 
We observe that this is indeed the case on our two finest 
lattices whereas for the coarsest, scaling violations are 
apparent. In the same figure we also show Am ob- 
tained using unquenched, dynamical Wilson configura- 
tions. It nicely falls on the same curve as the fine- 
lattice quenched results. This corroborates that scal- 
ing has set in and that quenching effects at these quark 
masses are small. Therefore a quenched study of diquark 
properties in the quark mass range used here is a very 
good approximation and will be adopted in the rest of 
the discussion. The Ansatz given in Eq. ^ provides 
a good fit to the lattice data. Using the lattice data 
with good scaling behavior at /3 = 6.0 and 6.2, we find 
C2 = 0.78 ±0.15 GeV, and ci = 0.138 ±0.01 GeV, which 
is the mass difference at the chiral limit [l^. This mass 
splitting can be compared to the A-nucleon mass split- 
ting, SruAN- Using our quenched data at /9 = 6.0 we find 
for the ratio Am/SniAN = 0.67(7), 0.73(8) and 0.67(8) 
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at K = 0.153, 0.154 and 0.155 respectively. Therefore our 
estimate of Am is consistent with the predicted 2/3 of 
the A-nucleon mass difference given in Ref . [8] . 
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FIG. 2: Mass difference between "good" and "bad" diquarks 
as a function of at /3 = 5.8 (crosses), (3 = 6.0 (open 
triangles), /3 = 6.2 (filled circles) and unquenched result (filled 
square). The dashed lines are fits to Eq. (|4}. 
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FIG. 3: Cr(r/a = 7.35,6*) versus cos(6l) for /3 = 6.2 and 
K. — 0.1520 for the "good" diquark ( asterisks) and the other 
diquark channels using the same notation is as in Fig. [T] Cor- 
relators for the "bad" diquark have been multiplied by two 
and for all negative parity channels by ten. 

Having determined the diquark spectrum, we now turn 
to the analysis of their structure by studying the density- 
density correlators defined in Eq. ([T|). The time t where 
the density operators are inserted is shown in Fig. [T] by 
the arrow (at t/a = 8), and is within the plateau range 
of the effective mass. We have verified that the correla- 
tors remain, within errors, unchanged when we vary the 
source-sink separation or equivalently t. Our aim is to 
fook for spatial correlations between the two light quarks 
in the various diquark channels. We take the location of 
the static quark as the reference coordinate from which 
the distances r„ and of the two quarks of different fla- 
vors are measured. We are interested in intrinsic diquark 
correlations, which persist as the static quark is moved 
away from the diquark. We therefore consider spherical 
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FIG. 4: Left: Cr(r/a = 5.1,e)/C-,^{r /a = 5.1,0) versus 
cos(6l). Right: Cr(r = 0.5 fm, r„d)/C.y5 (r = 0.5fm,0) ver- 
sus r„d, for the "good" (asterisks) and "bad" (filled triangles) 
diquarks at the lightest pion for our three lattice spacings. 

shells |r„| = \rd\ — r oi increasing radius r. Since the 
system is spherically symmetric the correlator depends 
only on r and the angle 9 — arccos(f„.fd) 0]. In the ab- 
sence of any correlation, the distribution of Cr{r,0) will 
be uniform as a function of cos(6'). Attraction will show 
up as an enhancement at small angles, near cos(0) = 1. 
Our cubic lattice breaks rotational symmetry and dis- 
torts the uniform spherical distribution, particularly at 
small angle 9. To remove such lattice artifacts, we nor- 
malize our distributions by a uniform lattice distribution. 
We show the resulting density correlators in Figs. [3] and 
|3]for various shell radii r. The physically relevant corre- 
lations are those that survive when r is large. We indeed 
observe that the qualitative behavior of the distributions 
does not depend on r, confirming that the color field gen- 
erated by the static quark does not affect the physics of 
diquarks once r is large enough. 

In Fig. [3] we show Cr{r, 9) for all the different channels 
as a function of cos(6') when the shell radius is fixed to 
r — 0.5 fm {r/a = 7.35). We clearly observe that the 
"good" diquark shows enhanced correlations at small 9, 
indicating attraction between the quarks. There is also 
a gradual increase in the distribution for the "bad" di- 
quark indicating a weaker attraction in this channel. This 
behavior persists as we increase or decrease the quark 
masses. For the negative parity channels, correlations are 
absent or very weak, except for the vector channel with 
r = 75 7i. However those disappear for lighter quark 
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masses. The negative-parity correlators are noisier be- 
cause the diquarks are heavier, so that a higher statistics 
analysis will be needed to fully determine the charac- 
teristics of their distributions. For this reason, we now 
restrict our analysis to the positive parity diquarks. 
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FIG. 5: String formation for three quarks, where one quark 
has been taken to the distant left. As a function of the sepa- 
ration between the other two quarks, the static potential rises 
with an effective string tension cJcB ~ o"-^- 

In Fig.|3]we show the density correlators for the "good" 
(scalar) and "bad" (vector) diquarks for all three lattice 
spacings at our lightest quark mass. The left hand set 
shows the correlator as a function of cos(0) normalized to 
1 at = for r / a = 5.1, which corresponds to r = 0.69, 
0.47 and 0.35 fm at /3 = 5.8, 6.0 and 6.2 respectively. 
As can be seen, irrespectively of the distance the "good" 
diquark shows stronger spatial correlations. Within this 
framework we can also extract the diquark size: At fixed 
r we look at the dependence of the correlator on the rel- 
ative u — d separation, r^d — 2r sin (0/2). In the right 
hand of Fig.[H we show the correlator (normalized to 1 at 
the origin) as a function of r^d for a fixed physical shell 
radius r = 0.5 fm, at our lightest quark mass. An ex- 
ponential dependence of C.y^{r,rud) oc exp(— rt((i/ro(r)) 
provides a gauge invariant definition of the diquark size 
ro(r) for a given value of r. The curves in Fig. [4l ob- 
tained from fits to an exponential dependence, describe 
well C^r^{r,rud) at all (3 values. Moreover the physical 
size ro(r) is the same on the two finer lattices, confirm- 
ing continuum- like behaviour as for the diquark masses. 
We can then evaluate rp as a function of the shell ra- 
dius r. Our measurements show a mild increase of rg 
with r, from ~ 0.9 fm to ^ 1.3 fm as r increases from 



0. 3 to 0.75 fm. For large r > 0.5 fm, rg is consistent 
with having reached a plateau. This suggests that the 
static source has a small influence on the "good" diquark 
giving a characteristic size of about 1.1 ± 0.2 fm. For 
comparison, at the same quark mass and using the same 
definition of size, we find for the p a size of 0.7 fm 9] . The 
large diquark size can be understood from the following 
qualitative argument: As the static quark is moved away 
from the diquark a q-q string tension develops. Consid- 
ering an increase e in the separation between the quarks 
in the diquark as illustrated in Fig. \E\ one finds that the 
effective q-q string tension is (-v/S/S) a where a is the q-q 
string tension. Since in the perturbative regime, the q-q 
attraction is also weaker than the q-q attraction, this time 
by a factor 1/2, the conclusion is that for all distances 
one expects the q-q attraction to be weaker. Thus a q- 
q diquark should be somewhat larger than a q-q meson, 
which is what we find. Size measurements for the "bad" 
diquark, on the other hand, show neither scaling nor con- 
vergence to a plateau value. The large values obtained, 
often similar to our box size, corroborate the weakness 
of spatial correlations in this channel. 

In conclusion, we have evaluated the mass splittings 
and density correlators of the complete set of diquark 
channels created by local diquark fields. Both observ- 
ables confirm the phenomenological expectation that 
QCD dynamics favors the formation of "good" diquarks, 

1. e. in the scalar positive parity channel. The character- 
istic size of this diquark, 0(1) fm, is large but consistent 
with the scale 0(200) MeV of the attraction. Even a 
"good" diquark is a large object, which may limit its rel- 
evance to hadron structure. The positive parity vector 
channel is higher in energy by about 2/3 the A-nucleon 
mass splitting, and forms an even larger object. All the 
negative parity channels have much higher energies. 
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